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ABSTRACT 

High-resolution observations of a quiet Sun internetwork region taken with the Solar 1-m 
Swedish Telescope in La Palma are analyzed. We determine the location of small-scale vortex 
motions in the solar photospheric region by computing the horizontal proper motions of small- 
scale structures on time series of images. These plasma convectively-driven swirl motions 
are associated to: (1) downdrafts (that have been commonly explained as corresponding to 
sites where the plasma is cooled down and hence returned to the interior below the visible 
photospheric level), and (2) horizontal velocity vectors converging into a central point. The 
sink cores are proved to be the final destination of passive floats tracing plasma flows towards 
the center of each vortex. We establish the occurrence of these events to be 1.4 x 1CT 3 and 1.6 



x 10 vortices Mm min respectively for two time series analyzed here. 
Key words: Sun: convection - Sun: granulation: - Sun: photosphere . 



1 INTRODUCTION 

The solar photospheric plasma is in constant evolution and is driven 
by convective processes in many spatial and temporal scales. Gran- 
ular convective motions are of great interest as this process of en- 
ergy exchange affects the evolution of magnetic structures while 
embedded in the plasma and might end up changing the topology 
of small emerging magnetic loops. By using high-resolution time 
series at rapid cadences we are nowadays able to register some of 
the finest and fastest scales of solar activity in an attempt to de- 
termine the Sun's structural building blocks. Investigations aim to 
describe the configuration and interactions of these features (i.e. 
magnetic concentrations) on their way up from the photosphere to 
upper layers in order to account for a more complete picture of solar 
activity. 

The quiet Sun has proved to be more dy- 
namic and magnetically active than prev io usly thought 
(Dommguez Cerdena, Kneer & Sanchez Almeida 2003), in- 
volving the rapid evolution of bright points (hereafter BPs) mostly 
associated to strong magnetic fields (up to the order of a few kG). 
Depending on their strength, magnetic concentrations are affected 
by the convective motions, being dragged and reconfigured 
through, for ins tance, fragmentation and coalescence phenomena 
(Viticchie 2009). Moreover, small-scale and short-lived emerging 
mag netic loops represent a substantial part of the total magnetic 
flux dMartmez Gonzalez & Bellot Rubiol 120091) and therefore are 
useful to understand how they interact with the photospheric 
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plasma. It is also important to remark that granular convection 
has been proposed to be involved in the mechanism driving the 
production of an effi cient turbulent dynamo re sponsible for the 
quiet Sun magnetism l |Petrovav~& Szakaly 1 9931) . 

Vortical motions in the quiet photosphere associated to in- 
tense downdrafts are commonly found i n theoretical simulations 
dSpruit. Nordlund & Title] ll99ok Izirkej 1 19931: IStein & Nordlund 
Il998t) . Very recent simulations by IShelvag et alj d201 ll) have also 
shown that large amount of vorticity in the photosphere is formed 
due to the interaction of plasma with the magnetic fields in the in- 
tergranular lane junctions. These types of vortex-like motions have 
been o b servationally repor te d at large spatial scal e s bv | Brandt et al.l 
dl988h : iTitle et al.l dl992h: lAttie. fanes & Pottsl 12009]) : at small 
sc ales by iBonet et alj 120081): and in the chromosphe ric quiet Sun 
bv lWedemever-Bohm & Rouppe van derVoortl(l2009h . 

In the present work we concentrate on a quiet Sun internet- 
work region that has been target of a few recent stud ies describ- 
ing the pre sence of convective dri ven vortex flows (Bonetetal.1 
120081 l20ld) . the coverage of BPs dSanchez Almeida et al.ll201(ih . 
and e vidence of magnetic c oncentrations dragged by a swirl mo- 
tion tealmaceda et alj|2010h . We use time series of G-band images 
recorded at high cadence in order to compute the proper motions 
of structures in the observed region. By generating the horizontal 
velocity maps we discriminate the regions displaying a pattern of 
converging flows into a central point of intense strong downdrafts 
where the plasma is returned to the solar interior. These events are 
then tracked within the field-of-view (hereafter FOV). We estimate 
their lifetime to be around 10 to 20 min. Our results evidence that 
the detected events can not be considered as homogeneously dis- 
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Figure 1. Example of one of the detected events displaying horizontal flows 
converging into a central point. Arrows represent the horizontal velocity 
vectors inferred from LCT. The displayed swirl motion exhibits a counter- 
clockwise sense of rotation. The background image represents the 20-min 
average image in false color with a FOV of 3"x 3". 

tributed over the FOV. We highlight the importance of understand- 
ing the effects of photospheric vortical motions on the configura- 
tion of magnetic fields in their way up to higher solar atmospheric 
layers, that are yet to be found. 



2 OBSERVATIONS AND DATA PREPARATION 

The data were acquired during a particular observing run on 29 
September 2007 wi th the Swedish 1-m Solar Telescope (SST, 
Schar mer et al .120031) in La Palma, Canary Islands, as part of a long 
international campaign with joint observations making use of other 
solar facilities at the Canary Islands Observatories. The region of 
interest corresponds to a quiet Sun area close to the solar disc cen- 
tre (p=0.99). Images in G-band (/1430.56 nm) were recorded at a 
fast cadence and the Multi-Frame Blind-Deconvolution (MFBD, 
Lofda hlll2002h restoration technique was applied to correct the im- 
ages from the aberrations induced by the turbulent atmospheric 
medium that affect their quality. The restored images with an ef- 
fective FOV of 69 x 69 square arcsec and spatial sampling of 
07034/pix were corrected with the standard procedures of flat- 
fielding, dark current substraction, elimination of spurious pixels 
and borders, and were grouped in two continuous time series of 
images, si: 08:47 to 09:07 UT, and s2: 09:14 to 09:46 UT, with 
15s cadence. The ~7-min gap between si and s2 corresponds to 
poor quality images that had to be discarded due to bad observing 
conditions. The final steps included the compensation for diurnal 
field rotation, destretching and subsonic filtering to eliminate resid- 
ual jittering. More deta ils on the data preparation can be found in 
IBalmaceda et aljfeoid) . 



3 SINKS AND CONVERGING FLOWS 

Our analysis is based on the widely used local correlation track- 
ing techniques (L CT, Nove mber & Simonlll988h . implemented by 
Molo wnv-Horas & Yil dl994h . We compute the horizontal proper 
motions of structures using a Gaussian tracking window with a 
Full- Width at Half-Maximum (FWHM) of l'.'O, to generate maps 
of horizontal velocity (flow maps). The same procedure was used 
by iBalmaceda et al. on the same data for the analysis of 



a fraction of the whole FOV (10" x 10") displaying a region of 
strong negative divergence. These authors found converging hori- 
zontal flows, i.e the velocity arrows pointing to a common destina- 
tion being the location of the sink. 

In this work we present the results obtained from the LCT 
analysis for series si and s2 over the complete FOV for 20-min 
intervals. From the horizontal velocity maps we compute the di- 
vergence field defined as: Vv = ^7 + ^7, where v. v and v y are the 
corresponding x and y components of the horizontal velocity vec- 
tor. 

Vertical velocities are inferred after multiplying the horizon- 
tal flow divergence by the so-called scale h eight of the flux mass 
(h m = 150 km), following iNovembed dl989h . Note that the term 
divergence refers to the divergence of the horizontal velocities. For 
a detailed derivation, physic al explanation and validity of this re- 
lation, we refer the reader to iMarquez. Sanchez Almeida & Bonen 
(2006). Resulting vertical velocities are manifestly conditioned by 
the LCT average and therefore can not be directly compared to 
Doppler velocities. An upcoming work attempts to establish mar- 
gins for the comparison of the two. In our data the velocity values 
as corresponding to a 20-min average range between 1.8 (upflows) 
and -1.2 km s _I (downflows). The overall distribution of horizontal 
velocities of one of these events is shown in Fig. [I] There is a clear 
trend in the direction of the velocity vectors, pointing towards the 
centre of the image (draining point) and forming a swirl motion 
with a counterclockwise sense of rotation. The background image 
is the G-band intensity average over the 20-min interval. The centre 
coincides with a dark structure formed by the junction of various 
intergranular lanes. 

Figure [2] displays the vertical velocity maps for si (upper 
panel) and s2 (lower panel) respectively, using a common time 
coverage of 20 min for both time series (i.e. the total duration of 
si and the first 20 min of s2) in order for the results to be compa- 
rable. Regions showing intense upflows denoted in different shades 
of red color correspond to intensive and recurrent activity from ex- 
ploding granules with velocities of ~1 km s , whereas strong sinks 
with negative, i.e. entering the plane of the figure, vertical veloci- 
ties of magnitude ~0.9 km s -1 are observed as dark features. The 
intergranular lanes are generally traced by the elongated structures 
in blue with downflows (-0.2 km s" 1 ). The black box i n the upper 
panel represents the FOV studied by IBalmaceda et alj < l20lbn . En- 
circled in black are the locations of the detected vortical motions 
(also referred to as events) from the horizontal flow pattern. The 
selection was performed by visual inspection, considering a vortex 
where velocity vectors in the flow map was converging and chang- 
ing direction defining a swirl. The area coverage of the swirl events 
are assumed to be circular and drawn as circles in the maps follow- 
ing the mentioned visual criteria. In all detected cases the events are 
located in regions with strong downflows (v ; — 0.8 km s ). Note 
that v. represents the averaged v LO s over the time elapsed. The pres- 
ence of all the events in at least two of the consecutive flow maps 
over 5-min intervals suggest that they remain coherent in a range 
between 10 and 20 min. The positions of the events change from 
the first to the second time series after the 7-min gap in the observ- 
ing run (dotted circles in the lower panel in Fig.|2]show the location 
of the events in the upper panel in the same figure for comparison.) 
These events do not seem to be homogeneously distributed over the 
FOV but rather grouped in certain locations resembling the meso- 
granular and supergranular patterns. 

Figure [3] shows the final destination of passive tracers 
or floats (white dots in the figure) driven by the computed 
mean field of horizontal velocities, as commonly used to trace 
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Figure 2. Maps of vertical velocities derived from the horizontal velocity field using LCT analysis for the SST/G-band ti me series si (upper pan el) and s2 
{lower panel). The same time coverage of 20 min is used in both time series. The black box represents the FOV analyzed by Balmaceda et al. (2010) including 
a pair of strong swirl events. Encircled in black are regions in which the horizontal velocities are converging towards a central point as inferred from the flow 
maps. The size of circles indicates the area coverage of the corresponding swirl event. Red circles denote peculiar regions with vertical velocities that do not 
conform vortical flows. Dashed circumferences in the lower panel show the location of the events in the upper panel for comparison. 
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Figure 3. The images show the final destination of passive tracers (in white) driven by the computed mean field of horizontal velocities. The background 
represents the high-contrast image smoothed by the 20-min average in every case. Locations of the swirl events detected in Fig. [2] (in black circles) are 
encircled in blue/yellow denoting clockwise/counterclockwise sense of rotation, respectively. Red circles are regions in which there is a strong concentration 
of floats, also plotted in Fig. [2] 
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the evolution of plasma motions (see for instance |^ 1 19921 ; 
Marquez, Sanch ezAlmeida & Bonetll2006l) . The background rep- 
resents the averaged image in every case. The position of the events 
detected in Fig. [2] are this time encircled in blue/yellow denoting 
clockwise/counterclockwise sense of rotation. For both time series 
we found a larger number of counterclockwise motions, though the 
difference is not very significant. In general these events are evi- 
denced by a larger concentration of floats with the exception of a 
very few cases in the second time series (lower panel) in which no 
sign of floats is detected though the areas are characterized by a 
junction of somewhat dark intergranular lanes. The distribution of 
vortices shows therefore a prominent correspondence with the lo- 
cation of intergranular lanes. The red circles in Fig. [3] are regions 
in which there is a strong concentration of float tracers, forming 
in some cases ring-shaped features. Although some of these re- 
gions display converging flows and correspond to downflows, they 
are not linked to very strong downflows when plotted over the im- 
age of vertical velocities (see red circles in Fig.[2j. Their location 
also change from the first to the second time series (e.g. compare 
the distribution of red and dotted-red circles in the lower panel of 
Fig. [2j. They do not represent though the most significant part of 
the detected events but ~10% in both series, respectively. We have 
calculated the occurrence of the detected vortex-type events in the 
FOV obtaining 1.46 x 1(T 2 and 1.09 x 1CT 2 vortices Mm" 2 for 
the SST/G-band series si and s2, respectively (the values increase 
up to 1.66 x 1CT 2 and 1.25 x 1(T 2 vortices MnT 2 if including the 
events not associated to very strong downflows and encircled in red 
in Fig.0) 



4 STATISTICS OF SWIRL MOTIONS 

As discussed in the previous section we have computed horizontal 
velocities and their corresponding divergences from G-band obser- 
vations. With the location of the detected vortices (circles in Fig. [3} 
we define a binary mask used to compute some statistical properties 
in vortex-type areas (i.e., a mask of vortices). The vorticity is a vec- 
torial magnitude whose direction is perpendicular to the XY plane. 
Using the equation: V x v = - we compute vorticities of the 
horizontal flow, with v- not included in the computation. Previous 
works have al ready studied vorticity in the solar p hotosphere from 
observational dWang at al.llT 995 : Bon et et al. l2010h and theoretical 
dStein & Nordlundlll998l : IShelvag et al.ll2011l) approaches. 

Figure [4] (upper panels) shows histograms for the whole FOV 
of ~ 50 x 50 square Mm. Panel a) displays independent results for 
series si (red) and s2 (blue) and total values for both series (black), 
resulting in a Maxwellian distribution of horizontal velocities, 
with a most probable value of 0.52 km s~'. Histograms in panel 
c) display a similar distribution of speeds using exclusively values 
within the mask of vortices. The most probable speed reduces to 
0.48 km s _1 . Median values combining si and s2, referred to as 
median(ji +J 2), are shown in each panel in Figure |4] Independent 
median values of horizontal speed for sl/s2 are 0.50/0.55 km s _1 
in Figure^ and 0.51/0.44 km s _I in Figure[4j:. Large horizontal 
velocities in the FOV are predominantly coming from recurrent 
exploding events (at the mesogranular scale) with large positive 
divergence values and hence the mean speed within vortex areas 
is expected to moderate as corresponding to negative divergences 
(convergences). 



Same as with horizontal velocities and divergences, we have 
also plotted vorticity values over the whole FOV and in the mask 



of vortices (i.e. in the areas inside circles in Figs. [2] and [3). Fig- 
ures^) and|4jl) show the distribution of divergence (dashed lines) 
and vorticity (solid lines). Independent results are again plotted for 
series si (red) and s2 (blue) and total values for both series (black). 
When considering the total FOV, the divergence distribution is cen- 
tered at a slightly negative value indicating a preference for con- 
verging flows probably located at the edges of supergranular cells. 
This behaviour is more evident when considering vortex areas only 
(panel d). The vorticity distribution for the whole FOV shows a 
slight trend for clockwise sense of rotation. Areas presenting vor- 
tices instead show a counterclockwise trend. 

Figure [5] plots a series of histograms computed using exclu- 
sively average values within the mask of vortices (i.e. mean values 
for each vortex). Vortices detected in series si and s2 respectively 
show different vertical velocities vz, with more intense downflows 
observed in the second one. From the histograms in Fig. |5^, most 
of the vortices are characterized by velocities of -0.4 km s and 
-0.6 km s _I for si and s2, respectively. It is important to bear in 
mind that vz is always negative since all vortices are observed in 
downflow areas, i.e. areas where the divergence is negative, with 
vz being proportional to divergence times the scale height. Fig. (5jj 
shows the distribution of the same quantities from Fig.|4^ and|4j; 
but considering the mean values of speed in areas displaying vor- 
tices. The results are similar. In Fig. [5]: and Fig.|5jl the distributions 
for positive and negative vorticities are shown. The medians for the 
counterclockwise and clockwise sense are comparable, indicating 
an apparent non-preferred sense of motion though in Section [3] we 
found counterclockwise events as being slightly more numerous. 
Fig. [5£ reveals a strong peak in the values between 240 and 250 
km for the mean radius. It is worth noting the small size of the vor- 
tices that remain stable during the two series since we would expect 
larger sizes i.e. of mesogranular scales, considering the 20-minute 
average used for the LCT computation. In Fig. [5j we include the 
circulation, defined for each vortex as the mean absolute vorticity 
times their corresponding area. A peak value of 1584.9 km 2 s -1 is 
clearly enhanced in the plots. 

Figure [6] displays scatter plots comparing different properties 
that characterize the vortices for series si (blue) and s2 (red). Lin- 
ear fits are shown in black for the combined values of both series. 
For all panels regressions were also separately computed for each 
series (not shown in the plots) with the resulting fits matching well 
between the two series. Panels a) and b) show the dependence of the 
vortex area and its corresponding velocities: larger vortices present 
higher horizontal and vertical speeds. This behavior is also clear 
in panel c). Larger scattering is observed in panel d), where the 
dependence of mean horizontal speeds with vorticities is shown. 
However, the linear fit shows a trend of slightly higher vorticities 
for larger speeds. 



5 DISCUSSION 

Recent advances in solar instrumentations from ground-based and 
satellite telescopes are providing us with unprecedented high- 
resolution data allowing the study of the dynamics of the finest 
solar structures. The evolution of the solar granulation pattern at 
the photospheric level is of great interest due to the interaction be- 
tween plasma and magnetic concentrations. In this scenario vor- 
tices, formed injunctions of multiple intergranular lanes, are places 
in which the concentration of these magnetic features are expected 
to be favoured an d therefore play an im portant role in the dynamics 
of the quiet Sun dKitiashvili et alj|201fj) . 
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Figure 4. Panels a) and c): Distribution of horizontal speeds. Panels b) and d): Histograms of the values of divergence (dashed lines) and vorticity (solid 
lines). Upper panels plot values in all the FOV whereas lower ones account for values within areas displaying vortices. All plots are computed using values 
for series si (in red), s2 (in blue) and combined values for both time series (in black). Corresponding most probable values for the combined case are shown 
in the panels. 



We search for the location of small-s cale swirls by trac king 
plasma motions in the same region where Bo net et al ](g00i) de- 
tected small whirpools by following the trajectories of BPs being 
swallowed by them. From 20-min average flow maps we identify 
events displaying a converging pattern of horizontal velocity vec- 
tors towards a central point that correspond to the strong sinks as 
initially predicted by the numerical simulations. These events are 
in all cases detected in the vertices of multiple granules along very 
intense intergranular dark lanes. This is clearly seen on the average 
image for the duration of the time series where the low intensity 
junctions are enhanced (in Fig. [3}. Data from two time series and 
common examination allow the comparison of results. A total of 70 
in si and 77 in s2 vortices have been detected, resulting in a density 
of 2.8 x 1CT 2 vortices Mm 2 and 3.1 x 1(T 2 vortices MnT 2 for si 
and s2 respectively. 

Averaging over the 20-min window we obtained space-time- 
density values of 1.4 x 10~ 3 vortices Mitt 2 min -1 and 1.6 x 
10~ 3 vortices Mirr 2 min~' for si an d s2 respectively. V alues are 
comparable to the number obtained by Bon et et ail d2008h of 1.8 x 
10~ 3 vortices Mirr 2 min 1 and lower than th e space-time-densit y 
of 3.1 x 10~ 3 vortices Mm" 2 mhr 1 found bv llonet et alj feOld) . 
It is necessary to bear in mind that the results from our flow fields 
are smoothed by the 20-min averages and the size of the employed 
tracking window. Many short-living vortical motions are likely di- 
luted by the LCT temporal average. 

BPs are distributed a l l ove r the FOV in this region as found 
bv lSanchez Almeida et al. I fcoioh from which some des cribe spiral 
trajectories whilst being engulfed by downdrafts that iBonet etHil 



d2008h described as convectively driven vortex flows. We have 
compared t he location of our detected vortices with the events dis- 
covered bv lBonet et al.ld2008l) in the same solar region. As a result, 
we have found 68% of coincidences but the remaining 32% corre- 
sponds to cases in the near vicinity of our vortices (i.e. less than 2" 
from the edge of the circular area enclosing the events in Fig.[2]and 

Rg. m 

We find for the series si and s2 counterclockwise/clockwise 
sense of rotation frequency value ratios of 53/47 and 52/48 %, re- 
spectively, which in view of the number of studied events shows 
no significant diff erence with e qually probabl e sense of rotation , 
in agreement with lBonet et alj ( 120081) . though IBonet et alj J201oh 
find a significant preference in favor of counterclockwise sense of 
rotation. The latter explain the different results assuming that the 
vortex rotation is influenced by the latitudinal so lar differential ro- 
tation, as the observations in iBonetet all [2010) were done at mid 
solar latitude, while the observations in lBonet et alj d2008l) and, of 
course, in this paper were done at the equator, where the influence 
of differential rotation is negligible. 

In te rms of dimensions o ur detected vortices coincide with the 
values in Bon et et al. Uioich of less than 500 km of radius. The 
majority of our detected vortices exhibit a radius of 241 ± 25 km. 
Some vortices, however, have a radius over 400 km that might not 
correspond to the main vortical scale detected but to signatures of 
larger scales of vortical motions. For the circulation we also found 
a pronounced peak value, at 1585 km 2 s . Thi s value is nea r ly half 
the 4000 km 2 s -1 for the vortex presented by iBrandt et ail dl988l) 
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Figure 5. Histograms for data within the mask of vortices showing mean values of: a) vertical velocity, b) horizontal speed, c) positive vorticity, d) negative 
vorticity, e) radius of vortices and f) circulation, computed for the vortices in series si (in red), s2 (in blue) and combined values for both time series (in black). 
Corresponding most probable values for si, s2 and for the combined case (sl+s2) are shown in the panels. 




Figure 6. Scatter plots of properties characterizing the detected vortices. Diamonds represent the values for si (red) and s2 (blue) . Black lines are the result 
of linear fits applied over the data with their corresponding regression coefficients and errors shown. 



which covers an area of about 5 granules and hence corresponds to 
a much larger mesogranular scale. 

Strong interacting vortices may play a role in excitation of 
solar acoustic oscillations, as shown bv lKitiashvili et all d20ld) . In 
particular these authors propose that high-speed vortices can attract 
and capture other vortices of opposite sense of rotation and their 
interaction might lead to their partial annihilation. In our sample of 
vortices we find several examples with these characteristics, though 



a deep analysis of the interaction of opposite-sign vorticities is out 
of the scope of the present work. 

More observational evidence from spectropolarimeric data 
would be required to complete the analysis and determine up 
to what extent the magnetic concentrations in the nearest vicin- 
ity of the detected swirl events are affected by the appearance 
of the vortex. An example of magnetic features appearing as be- 
ing dragged by this typ e of convective motions was presented by 
Balma ceda et alj fcOld) in one of the vortices within the FOV. This 
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type of vortical convective flows are thought to contribute gath- 
ering magnetic fragments thus amplifying weak magnetic fields 
near the surface and representing an important mechanism for the 
forma tion of more stable features dDanilovic. Schiissler. & Solankil 
l20ld) . Exhaustive analysis of the intensification of magnetic field 
in vortex-like regions and its morphology at different solar atmo- 
spheric layers will also contribute with precise inputs to improve 
numerical simulations development. 
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